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The static structure factors of the XXZ model in the presence of uniform field are determined 
from an exact computation of the groundstates at given total spin on rings with A'^ = 4, 6, . . . , 28 
sites. In contrast to the naive expectation a weak uniform field strengthens the antiferromagnetic 
order in the transverse structure factor for the isotropic case. 
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10 ■ I. INTRODUCTION 

o 

The antiferromagnetic properties of the one-dimensional spin 1/2 XXZ-model with Hamiltonian: 

H = 2Y,[Siix)Si{x + l) + S2{x)S2{x + l)+cosj S3{x)S3{x + l)] (1.1) 

x=l 

have been studied by analyticalE and numericali^^ methods. The critical exponents r]j,j = 1,3, which govern the 
large distance behavior of the spin-spin correlators in the groundstate: 

'gl {0\S,{0)S,ix)\0) ^ i = l,3. 



are given by:E 



X ■ 771 = ??3~i = 1 - ^ for < 7 < TT. 

u ■ 



At finite temperature T the spin-spin correlators decrease exponentially with a rate given by the inverse correlation 
length :eI 

^'(7,^)^$(7).T^ 

with a critical exponent: 

1^ = 1, 

independent of 7 and: 

sm7 V TT/ 

It has been shown in Ref. ^ that the structure factors: 

N/2-1 

Sj{-f,p=2TTk/N,T,N) = 1 + (-l)'^'4(0|5j(0)S'j(A^/2)|0) +8 ^ {0\Sj{0)Sj{x)\0) ■ cos(px). 
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are most suited to extract the critical behavior from finite systems. In Refs. 
factors in the following three limits: 

p = 7r, T = 0, N^oo, 



we have studied the static structure 



P ^ TT, T = 0, N = oo, 



p — TT, r — *■ 0, N = oo. 
The transverse and longitudinal structure factors have a common form in each of these limits: 

where: 

N 



(1.2) 



VN 



Vp = — — ' VT 



are the "running variables" . The common form reflects the fact that the structure factors scale in the critical regime: 

TT, T — >0, N — >oo, (1.3) 



if we keep fixed: 



P 



Zl 



N, Z2 



N 



C(7,r)- 



(1.4) 



The antiferromagnetic properties of the model become visible in the singularities of the structure factors in the limit 
(1.3). The transverse structure factors arc infinite in this limit; they develop a "hard" singularity. In contrast, the 
longitudinal structure factors stay finite. T heir critical behavior is hidden in subleading terms which produce a "soft" 
singularity. Away from the critical regime (O) the antiferromagnetic properties are lost. 

There are further possibilities to destroy antiferromagnetic ordering - e.g. by a uniform external field h or by 
frustration - i.e. by switching on a next-to-nearest neighbor interaction. In this paper we are going to study the 
effect of a uniform magnetic field on the zero temperature structure fa ctor s. For this purpose we have determined the 
static structure factors in the groundstates {Ss) of the Hamiltonian (1_T) at given total spin S3. The groundstates 
were computed with a Lanczos algorithm on rings with = 4, 6, 28 sites. __ 

We will discuss the static structure factors as function of the magnetization M — S3/N. The knowroLl magnetization 
curve translates M into the fieldstrength h. For convenience, we list in Table I the /i-field values corresponding to our 
M-values: 

TABLE I. magnetization and the corresponding ft-fields 



M 


1 


1 


1 


1 


3 


1 


« 


fi 


4 




s 




h{j = 0) 


0.96 


1.20 


1.59 


1.83 


1.91 


2.00 


/i(7 = O.Itt) 


0.93 


1.17 


1.54 


1.78 


1.86 


1.95 


/i(7 = 0.27r) 


0.83 


1.06 


1.41 


1.64 


1.72 


1.81 


/i(7 = O.Stt) 


0.38 


0.50 


0.71 


0.87 


0.92 


1.00 



The outline of the paper is as follows. In Sections II and III we discuss the characteristic features of the longitudinal 
structure factors as there are the p and M dependence and the finite-size effects. The same is done for the transverse 
structure factors in sections IV and V. 
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II. THE LONGITUDINAL STRUCTURE FACTORS AT FIXED MAGNETIZATION 



The lougitudinal structure factor of the XX-model (7 — tt/2) is known from the exact solution obtained by 
NiemeijerQ in the fcrniion representation: 

^ f 10 AT M\ 2 for 0<P<Pz{M) 

S3 7 = 7r/2,p,M,iV^ - <^ (2.1) 

L Ps^M) for P3{M) <p<n, 



where 



P3(M) = 7r(l-2M). (2.2) 



Though the resuh of Ref. || has been derived for the therniodynamical hmit iV ^ cjo, ( |2.lD turns out to be correct for 
all system sizes with = 4, 6, 8.... The linear behavior in p has been found beforecl for the case M = 0. 

Let us next turn to the isotropic case 7 = 0. Fig. 1 presents a panoramic view on the longitudinal structure factor 
with N — 20, 22, 28. The emergence of the singularity at p = tt, M = is clearly visible. Along the momentum 
axis at M = we see the logarithmic singularity: 

83(7 = 0,p,M = 0) = r-s{0) • In (1 - I) , (2.3) 

discussed in Ref. ^. Along the magnetization axis: 

p = Po = n, M = S3/N 0, (2.4) 
we observe a logarithmic singularity in M . The same type of singularity is also found in the limit: 

p = P3(M), M^O, (2.5) 



where p3{M) is given in ( |2.2| ). The longitudinal structure factor has its maximum aX p = p^{M)^ M fixed. With 
increasing strength of the uniform field the maximum position moves frorap = tt to p = - i.e. from antiferromagnetic 
to ferromagnetic order. Such a-behavior has been conjectitiKed by MiillerQ et al.. Indications for this have beeu-found 
also by Parkinson and BonnerE2l and Ishimura and Shiballil on small systems [N < 14). Johnson and FowlcrEJ were 
able to reformulate the isotropic Heisenberg model for large spins and magnetizations close to saturation in terms of 
a gas of magnons. By accident their prediction for the longitudinal structure factor in the limit M ^ 1/2 is identical 



with the exact result (2T) for the XX-model. In Fig. 2 we compare the longitudinal structure factors for 7 = and 
7 = 7r/2 at M = 1/4 and M = 1/3. For M = 1/3 the structure factors almost c oinc ide. For smaller M values, 
however, the isotropic structure factor deviates from ( |2.l| ). The cusp along the line ( |2.2| ) becomes more pronounced 
with decreasing values of M. 

Looking at the finite-size effects which will be analyzed in the next section we find an N~'^ behavior away from the 
cusp and a slower decrease N^^^ with ^3 « 0.5 at the cusp p = p3(M). The change in the finite-size behavior signals 
the emergence of a nonanalytic behavior in the thermodynamical limit. 

We have also determined the longitudinal structure factors for the anisotropics j/tt = 0.1, 0.2. The p—M dependence 
of the longitudinal structure factor looks similar to the isotropic case. Instead of an inf inity we find a peak at 
p — n, M = 0. In the limit p — > tt, M = the structure factor is adequately described by ( |l.2| ) with a — p. In the 
limit M Q,p = n we find again a behavior of the form (|1.2|) with a running variable: 



The appearence of the cusp along the line p — p3(M) is indeed independent of the anisotropy parameter 7. For 
increasing values of 7 the cusp is less pronounced. 



III. FINITE SIZE ANALYSIS OF THE LONGITUDINAL STRUCTURE FACTOR 



In the critical regime 



N^^, M-^^0, T^O, (3.1) 
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we expect the longitudinal structure factors to obey finite-size scaling: 

S3(7,P, M = S3/N, T, N) - 33(7; ^1, ^2, Z3) ■ S3(7,P, M, T, iV = ex.) 
In the combined limit ( |3.l[ ) we have to keep fixed zi, Z2 - defined in ( |l.4| )- and: 

Z3 = M-N^ S3. 
In Ref. ^ we checked finite-size scaling in the combined limit: 

T — > 0, N — > 00, Z2 fixed and p = tt, M = 0. 



(3.2) 



It was found that finite-size scaling works for the transverse structure factors at all 7 values. In contrast, finite-size 
scaling breaks down for the longitudinal structure factors, if 7 > O.Stt 

In this section we are going to study consequences of the ansatz 
The behavior of the isotropic structure factor can be read of (\ 



S3(7 = 0,p = ^,M,iV = oo) 



) in the limits (2.4) and ( p.SD at zero-temperature. 
M 



f 3(0) -In 



where i = 0, 3 stands for the two limits (2.4 275| ). The isotropic structure factors on finite systems with N sites scale 
in M (for M 7^ 0) and are linear in — InAf in the limit M 0, as can be seen from Fig. 3. Finite-size effects are 
small for po — n but rather large for ps = Tr(l — 2M). We have studied the finite-size effects in the differences: 



A, (7, M, N) = 83(7 = 0,p„ Af, N) + ln(2M), j = 0, 3, 



(3.3) 



at fixed magnetizations: 



M = 



1 3 

8' 8 

1 

4 

1 2 

6' 6 



N = 8,16,24 

N = 4,8,12,16,20,24,28 

N = 6,12,18,24, 



which can be realized on the systems with size N. The A'^-dependence of the difference (|3.3|) can be parametrized by 



Ao(7, M, N) = Ao(7, M) + co(7, M) ■ N'^, 
A3 (7, M, N) = A3 (7, M) + C3(7, M) ■ N-'^ 



(3.4) 



0.5. 



In our finite-size analysis we have included as well the magnetizations: 



M 



J_ _2_ J_ J_ 

10' 10' 10' 10 
J_ _5_ 
12' 12 

J_ J_ _3_ J_ _5_ _6_ 

14' 14' 14' 14' 14' 14 



N = 10,20 
N = 12,24 
N = 14,28, 



which occur on two systems. Here we have assumed that the finite-size dependence is described correctly by ([3.4|). 
The extrapolations of the structure factors to the thermodynamical limit are represented by the the solid dots in Fig. 
4. 

So far our estimates of the thermodynamical limit are restricted to magnetizations M > 1/14 due to the finiteness 
of our systems N < 28. Additional information on the structure factors for Af-values: 



M = 



1 

iV' 



N = 14, 16,. ..,28, 



(3.5) 



closer to the critical point M — can be obtained if finite-size scaling (3^) holds for the lon gitu dinal structure factors 
in the limits (p^.2.5) where we keep fixed 23 = 6*3. This variable is just 1 for the sequence (3.5), and we can compute 
the scaling function at this value from: 



ffi(7,Z3) 



^ S3(7^0,P = Pj,M^ l/14,jV= 14) 
^3=1 S3(7 = 0,p = Pj,M = l/14,7V = oo)' 



In this way we get from the finite-size scaling ansatz (3.1) in the limits (2.4,2.?;) an estimate of the thermodynamical 
limit of the structure factors: 
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S3(7 = 0,p = p„Af,iV = «3) = 



S3{j^0,p^p,,M^l/N,N) 
5j(7,l) 



J =0,3, 



(3.6) 



for the sequence of M values in ( |3.5| ). The result for the differences (3^) is marked by the open dots in Figs. 4(a), (b). 

We have repeated the finite-size analysis - described above for the isotropic case - for: 'y/w = 0.1, 0.2. The resulting 
estimate of the thermodynamical limit: 



versus the variable: 



Aj(7, M, TV = (x.) = S3(7,p„ Af, iV - (x.) - ^3(7, Af), j = 0, 3, 



is represented in Figs. 4(a), (b) by the triangles and squares, respectively. 



IV. THE TRANSVERSE STRUCTURE FACTORS AT FIXED MAGNETIZATION 



The most remarkable property of the transverse structure factor is its approximate constancy: 



Si(7,p,A//,r = 0,7V) « 2A/, for 0<p<2TrM, 



(4.1) 



which has been found by Miilleia et al. on small systems for the isotropic case 7 = 0. For S'3 = MN = N/2 — 1 
(4.1) can be e asily proven to be exact by means of the Bethe ansatz solution. For S'3 < N/2 — 1 and < p < 2Af7r, 
however, ( |4.l| ) is not exact. As an example we present in Figs. 5(a), (b) the momentum dependence of Si at Af = 
1/4, 1/3, N = 8, 12, ...,28 for 7 = 7r/2 and 7 = 0, respectively. The constancy in the regime p < 2ttM is striking. 
Deviations from (4.1) can be seen on a scale magnified by a factor 100 in the inset of Figs. 5(a),(b). These deviations 
follow a single scaling curve which increases monotonically with momentum p. On the scaling curve finite-size effects 
die out with N^^. At p = 2Af7r, however, significant finite-size effects of the order N^^^ with Si ~ 1.0 for 7 = and 
(5i ^ 1.3 for 7 — 7r/2 b ecom e apparent. 

Beyond the regime (4.1) the transverse structure factor of the XX-model (7 = n/2) is linear in (1 — p/vr)^^/^, as 
can be seen from Fig. 5(a). Therefore we find the same type of singularity for p ^ tt at Af = 0, 1/4, 1/3. In contrast 
to the naive expectation antiferromagnetic order is not destroyed in the transverse structure factors by an external 
field. 

In the isotropic case (7 = 0) the transverse structure factor is approximately linear in — ln(l — p/tt) for p > 2Af7r, 
as can be seen from Fig. 5(b). This type of singularity was found for p — > tt at Af = 0. The slight curvature at 
M = 1/4 might indicate that the type of the singularity has changed here to a power behavior (1 — p/7 r)~°. 

In the anisotropic case with 'y/n — 0.1,0.2 and M = 1/4 we find agai n the constant behavior (4.1) and a linear 
increase in (1 — p/7r)''i^^. Again this type of singularity follows from (1.2) with a—p and M = 0. 



V. THE TRANSVERSE STRUCTURE FACTOR AT CRITICAL MOMENTUM 

Let us start with the XX-model (7 — 7r/2). Fig. 6(a) shows the transverse structure factor for p = tt and 
N = 8, 10 , 28 as function of M. In contrast to the longitudinal case we have no scaling in M. At Af = we know 
from ( |l.2| ) with a — N that the transverse structure factors diverge as ^/N. The same type of divergence appears at 
M = 1/4, as can be seen from Fig. 7(a). Here we compare the A^-dependence of 81(7 — 7r/2,p = tt, M, N) for M = 
and Af = 1/4. In both cases there is a linear increase in y/N. 

Let us now turn to the isotropic case (7 = 0). Here the longitudinal and transverse structure factors coincide 
provided that there is no external field. In the presence of a uniform field, however, they differ drastically. At p = tt, 
the longitudinal structure factor is infinite at Af = 0, but becomes finite and monotonically decreasing for M > 0. 
In contrast the transverse structure factor stays at infinity for Af > 0. Moreover, on finite systems, the longitudinal 
structure factor scales for Af > - as was demonstrated in Fig. 3 - whereas the corresponding transverse structure 
factor (at p = tt) does not scale at all, as can be seen from Fig. 6(b). More surprising, at fixed Af > not too 
large, the transverse structure factors increase with the system size N stronger than at Af = 0. A comparison of the 
A^-dependence at Af = and Af = 1/4 is shown in Fig. 7(b). At Af =1/4 the increase with In is definitely steeper, 
which signals a strengthening of the singularity at p = tt. Note also that there are deviations from linearity in InA'^, 
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which increase with M . This might indicate a change from a logarithmic behavior at Af = to a power behavior for 
M > 0. 

In Ref. ^ Miiller et al. reported on the transverse structure factor for N = 10; they found already that the dominant 
mode remains situated at p — n independent of the field. Fig. 6(b) tells us that the strengthening of the singularity 
at p = n becomes more and more pronounced with increasing system size TV. 



VI. CONCLUSIONS 



In the presence of a uniform external field in z-direction the static structure factors of the XXZ-model show up the 
following features: 



1. The longitudinal structure factors have a cusp along the line (2_^). In case of the XX-model (7 ~ 7i'/2) there are 
no finite-size effects and the longitudinal structure factor is given by (2.1) for all system sizes with iV = 4, 6, .... 
For smaller values of 7 and M, the cusp becomes sharper. Finite-size effects decrease along the cusp with N~^^, 



63 ~ 0.5. Away from the cusp we find a more rapid decrease of the order N' 



-2 



The longitudinal structure factor is finite for 7 7^ and p ^ tt, M ^ 0, but develops a logarithmic singularity 
in this limit for the isotropic case (7 = 0). This means: A uniform field weakens the antiferromagnetic order 
in the longitudinal structure factor for 7 = 0. 



3. The transverse structure factor is almost contant for p < 2A/7r [cf. ([4.l| )]. Finite-size effects die out slowly with 
N~^^ with Si ~ 1 along the line p — 2M'k, but rapidly with away from this line. 

4. In the limit p ^ n, M =1/4 fixed, we observe a singularity of the type (1 — p/ir)'^^"^ in the transverse structure 
factor for 7 = 7r/2. In the isotropic case (7 = 0) this singularity appears to be stronger than — ln(l — p/tt). 
This means: A weak uniform field strengthens the antiferromagnetic order in the transverse strucure factor 
for 7 = 0. 

Therefore the effect of a uniform external field on the longitudinal and transverse structure factors for (7 = ) is 
similar to the effect of switching on the anisotropy parameter 7. The logarithmic singularity found in the isotropic 
structure factor at p = n changes with 7. It is strengthened in the transverse, but weakened in the longitudinal 
structure factor. 
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FIGURE CAPTIONS 

FIG. 1 The longitudinal structure factor versus momentum p and magnetization M with a ridge along the line p^i^M) = 
7r(l - 2M), for N = 20, 22, ...,28. 

FIG. 2 Comparison of the longitudinal structure factors at 7 = 7r/2 [Eq. and 7 = for M ~ 1/4 (o),l/3 (•), 
respectively. 

FIG. 3 Scahng of the longitudinal structure factors 83(7 — 0,p, M, N), j = 0, 3 with po — tt (open symbols), ps = 
7r(l — 2M) (solid symbols), with the magnetization M. 



FIG. 4(a) Estimate of the thermodynamical limit for the difference Ao(7 , M, N = 00) [Eq. (3_^)] at momentum po = tt: 
Solid symbols represent results from the finite-size analysis (3.4), open symbols are results from finite-size scaling 



FIG. 4(b) Same as Fig. 4(a) for the difference A3 (7, M, TV = 00) at the cusp p3{M) = 7r(l - 2M). 

FIG. 5(a) The momentum dependence of the transverse structure factor of the XX model (7 — 7r/2) at fixed magnetizations 
M =1/4 (o), M = 1/3 (•), respectively. The inset shows a magnification of the low momentum regime. 

FIG. 5(b) Same as Fig. 5(a) for the transverse structure factor in the isotropic case (7 = 0). The inset shows a magnification 
of the low momentum regime. 

FIG. 6(a) The transverse structure factor of the XX model (7 = 7i'/2) versus magnetization M at momentum p = tt. 

FIG. 6(b) Same as Fig. 6(a) for the isotropic model (7 = 0). 

FIG. 7(a) Comparison of the size dependence of the transverse structure factors of the XX model (7 = tt/2) aX p — tt: 
M = 1/4 (o),7\f = 1/3 (•), respectively. 

FIG. 7(b) Same as Fig. 7(a) for the isotropic model (7 = 0). 
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